arXiv:l506.00936v2 [hep-th] 17 Jul 2015 


Cosmological Attractors and Initial Conditions for Inflation 


John Joseph M. Carrasco, Renata Kallosh and Andrei Linde 
Department of Physics and SITP, Stanford University, 

Stanford, California 94305 USA 

Inflationary a-attractor models in supergravity, which provide excellent fits to the latest obser¬ 
vational data, are based on the Poincare disk hyperbolic geometry. We refine these models by 
constructing Kahler potentials with built-in inflaton shift symmetry and by making a canonical 
choice of the goldstino Kahler potential. The refined models are stable with respect to all scalar 
fields at all a, no additional stabilization terms are required. The scalar potential V has a nearly 
Minkowski minimum at small values of the inflaton field tp, and an infinitely long dS valley of 
constant depth and width at large tp. Because of the infinite length of this shift-symmetric valley, 
the initial value of the inflaton field at the Planck density is expected to be extremely large. We 
show that the inflaton field tp does not change much until all fields lose their energy and fall to the 
bottom of the dS valley at large tp. This provides natural initial conditions for inflation driven by 
the inflaton field slowly rolling along the dS valley towards the minimum of the potential at small 
ip. A detailed description of this process is given for a-attractors in supergravity, but we believe 
that our general conclusions concerning naturalness of initial conditions for inflation are valid for a 
broad class of inflationary models with sufficiently flat potentials. 


1. INTRODUCTION 


During the last two years, a new class of inflationary the¬ 
ories have been discovered: “cosmological attractors.” This 
class is very broad, including conformal attractors 0 , uni¬ 
versal attractors with non-minimal coupling to gravity |2], 
and a-attractors 00, and it also incorporates many pre¬ 
viously existing models such as Starobinsky model [7], GL 
model and Higgs inflation model m- Despite very 

different origins, all of these models make very similar cos¬ 
mological predictions providing an excellent match to the 
latest cosmological data mm- Moreover, these models 
can be further extended to describe not only inflation, but 
also the theory of dark energy/cosmological constant and 
supersymmetry breaking ISI- 

On a purely phenomenological level, the main features of 
inflation in all of these models can be represented in terms 
of a single-field toy model with the Lagrangian mug 


V=g 2 2 (i_ 


6ct 


- VM) . 


(i.i) 


Here 4>{x) is the scalar field, the inflaton. The parameter a 
can take any positive value. In the limit a —i oo this model 
coincides with the standard chaotic inflation models with 
a canonically normalized field <f> and the inflaton potential 
V{ct>) pH- 

However, the field <fi in © is not canonically normalized. 
It must satisfy the condition </> 2 < 6a, for the sign of the 
inflaton kinetic term to remain positive. One can easily go 
to canonically normalized variables tp by solving the equation 
= dip, which yields 

1— 


In terms of these variables, a tiny vicinity of the singular 
boundary of the moduli space at 4> = VGa stretches and 
extends to infinitely large ip. As a result, generic potentials 
V((f>) = V(\/6a tanhA|=) a t large tp approach an infinitely 
long dS inflationary plateau with the height corresponding 
to the value of V (</>) at the boundary: 

V dS = V(<j)) | 0=±v ^ . (1.3) 

This universal origin of the shape of these potentials leads 
to universality in their predictions, as explained in 000 
and formulated in a particularly general way in [14] : The ki¬ 
netic term in this class of models has a pole at the boundary 
of the moduli space. If inflation occurs in a vicinity of such 
a pole, and the potential near the pole can be well repre¬ 
sented by its value and its first derivative near the pole, all 
other details of the potential far away from the pole (from 
the boundary of the moduli space) become unimportant for 
making cosmological predictions. In particular, the spec¬ 
tral index depends solely on the order of the pole, and the 
tensor-to-scalar ratio relies on the residue m- All the rest is 
practically irrelevant, as long as the field after inflation falls 
into a stable minimum of the potential, with a tiny value of 
the vacuum energy, and stays there. 

This new class of models accomplishes for inflationary the¬ 
ory something similar to what inflation does for cosmology. 
Inflation stretches the universe making it flat and homoge¬ 
neous, and the structure of the observable part of the uni¬ 
verse becomes very stable with respect to the choice of initial 
conditions in the early universe. Similarly, stretching of the 
moduli space near its boundary upon transition to canonical 
variables makes inflationary potentials very flat and results 
in predictions which are very stable with respect to the choice 
of the inflaton potential. 


0 


= \/6a tanh 



( 1 . 2 ) 


The present work will pursue two different goals simulta¬ 
neously. First of all, we will study geometrical properties 
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of the moduli space in the supergravity realizations of these 
models, following nanzi- We will reformulate these models 
in terms of Kahler potentials and field variables which keep 
their geometric properties manifest. This new formulation 
will allow us to approach the problem of initial conditions for 
inflation in these models in a novel, more transparent way. 


The problem of initial conditions in these models is not 
quite trivial. In the simplest chaotic inflation models such 
as inflation may start very close to the Planck density. 

According to [16] [18ti22j . this makes initial conditions for 
inflation quite natural. However, in the new class of mod¬ 
els discussed above, as well as in the Starobinsky model and 
Higgs inflation, the inflationary regime begins at the energy 
density 10 orders below the Planck energy density. A solu¬ 
tion of the problem of initial conditions in such models was 
discussed in [23] . Here we will revisit it; we will show how 
this problem can be solved in the supergravity realizations 
of a-attractors. Most of our conclusions will have more gen¬ 
eral validity, being applicable, in particular, to generic non- 
supersymmetric attractor models (1.1). We will also show 
that in some cases, such as supergravity a-attractors with 
a C 1, inflation can begin at the density approaching the 
Planck density, thus reducing the problem of initial condi¬ 
tions to the one already addressed in [I61[I8H22]. 


There are two types of technical improvements of our a- 
models which we will develop in this paper. The first one, 
following mm, will allow us to use the Kahler frame where 
the inflaton shift symmetry is present in the new Kahler 
potentials. The second improvement with respect to earlier 
models corresponds to changing the Kahler potential for the 
goldstino multiplet, making it canonical rather than part 
of the logarithmic structure, which has a consequence: an 
improved manifest stability. 


We will make a choice of variables in which the infla¬ 
ton forms a Killing direction of the moduli space geometry. 
Namely, our holomorphic disk variable Z and the half-plane 
variable T used in HI 120 will be represented by the Killing 
adapted moduli space coordinates 


Z = 


T - 1 
T+ 1 


= tanh 


p + id 

V6a 


(1.4) 


Here the inflaton ip and the orthogonal field $ form a geom¬ 
etry independent on a Killing direction p: 


g w {^) = gmfd) = 



(1.5) 


As a result, the time evolution in our models with initial large 
kinetic energy, when the role of the potential is negligible, 
will be controlled by the fact that the momentum in the 
inflaton direction is preserved, namely 


P v = 0 where P v = a 3 {t)g w {f})p (1.6) 


This geometric fact helps us to argue that the total shift of 
the field p due to its initial velocity is about 10 Planck units 


or less, after which all memory about the initial velocity of 
the field tp at the Planckian time completely disappears. 

We will also numerically solve the Friedmann equations in 
FRW space-time metric for generic initial values of i9, p, p 
confirming our analytic analysis: we have an inflationary 
attractor behavior, where the memory about initial values 
of i9, p, p disappears and period of slow-roll inflation at 
the minimum of the potential at i9 = 0 takes place. 

We will show that with new Kahler potentials which have 
the inflaton shift symmetry in Z or T variables, the superpo¬ 
tentials are simpler and the relation between models in disk 
and half-plane variables simplifies. 

We will than proceed with the analysis of initial condi¬ 
tions for inflation in these models, with our new choice of 
variables, by making choices of initial values of the inflaton 
and its partner and by studying the time evolution of these 
fields, before and during inflation. The geometric nature of 
our models, and the existence of infinite dS valleys of con¬ 
stant width in our potentials, help to resolve this issue and 
allowing us to argue that the vast majority of initial condi¬ 
tions in these models leads to successful inflation. 


2. FROM DISK TO HALF-PLANE: NEW KAHLER 
POTENTIALS 


The cosmological attractor models can be described either 
in disk or in half-plane variables urn- The corresponding 
boundary of the moduli space, which plays an important role 
in these models, is either at ZZ < 1, or a half-plane with 
T + T > 0. 


Here we summarize the relation between disk and half¬ 
plane variables for generic case of 2-superfield models with 
our choice of the Kahler potentials and most general super¬ 
potentials. 


The relation between the Kahler potentials and superpo¬ 
tentials in the disk and half-plane variables requires a simple 
Caley transform, as suggested in [5j 


-—- 

T+ 1 


1 + 

1 - Z ' 


( 2 . 1 ) 


We will represent the Kahler potential in the following form: 


Kjs 


3a T (1 - ZZ f 

~2 ° S [{1- Z 2 )(l- Z 2 ) 


+ SS, 


( 2 . 2 ) 


W a = A(Z) + SB(Z) . (2.3) 

where S is a supermultiplet with a goldstino fermion and a 
sgoldstino scalar. This field may either belong to the usual 
unconstrained chiral multiplet, or it may be a nilpotent su¬ 
perfield as studied in m We will discuss both options in 
this paper. 
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We present our models as functions of a complex variable 
$ = ip + id , (2.4) 


where ip will be the inflaton in cosmological models and d 
will describe the orthogonal direction and 


Z = t anli —z —* (2.5) 

\/6a 


Our Kahler potential (2.21 in these variables has a manifest 
inflaton shift symmetry, ip' = ip + c 


/4> - $\ 
V \!/ 


+ SS . 
( 2 . 6 ) 


Ko(Z,Z) => K = —3a log cosh 
The superpotential is now 

Wb => W = A( tanh -=) +Sb( tanh -=) . (2.7) 


Note that in our models d = 0 during inflation and therefore 
the new holomorphic variable 4? during inflation becomes a 
real canonical variable. This is also easy to see from the 
kinetic terms in these variables, which are conformal to flat, 


ds 2 = 


dip‘ 


d-d 2 1 , 2 1 


— nQipipdf + n ppH)dd . 
2 / 2 „a l Z 


( 2 . 8 ) 


2 cos 2 \ if-d 

w 3a 


At i? = 0 they are canonical 


ds 2 


#=o — 


dip 2 + dd 2 


(2.9) 


We can also use the half-plane variables T + T > 0 where we 
have 


Km 



(T + T) 2 
ATT 


+ SS, 


( 2 . 10 ) 


= G(T) + SF{T) . (2.11) 


Now the disk and the half-plane models are related simply 
by the Caley transform (2.1), so that transition from one 
picture to the other is a simple substitution 


/ T — 1 - T — 1 \ . —. , 

Ko ( Z = T^l' Z = T^l) =i ^( T ’ T ) ■ < 2 - 12 ) 

and 

Wo( z = = ^h(T) . (2.13) 

This also means that 

G(T) = A(Z , F(T) = B(Z^). 

(2.14) 

When we hold SS outside of the log part of the Kahler po¬ 
tential, the field S does not change from one picture to the 
other. However, for any models with SS inside the log part 


of the Kahler the potential which we used before, the rela¬ 
tion between the goldstino multiplets in Z and T variables 
involves the dependence on the inflaton superfield, as shown 
in ;5j. We will explain below that when the field S is outside 
the log in the Kahler potential the inflaton partner is sta¬ 
ble for any a. Therefore we will focus here on models with 
canonical Kahler potentials for the S field as shown in eqs. 
(2221), (E6l) and pTlOb. 


3. a-ATTRACTORS AND THEIR STABILITY 


We will begin with a rather simple and general class of a- 


attractors in disk variables, with the Kahler potential (2.2) 
and superpotential 


w B = V^d sf(z). 


(3.1) 


Investigation of this theories simplifies considerably if during 
and after inflation the field S vanishes, along with the imag¬ 
inary part of the field Z. Indeed, as explained in El, the 


Kahler potential (2.2) has a shift symmetry under the shift 


of the inflaton field during inflation, when x = Im Z = 0: 
The Kahler potential vanishes independently of the value of 
the inflaton field z = Re Z. 


In that case, one can show that the potential of the canon¬ 
ically normalized inflaton field ip, which is defined by the 
relation 2 = tanh , is given by 


V = ap 2 f 2 (z) = ap 2 f 2 (tanh —)== 

V6a 


(3.2) 


The potential has an infinitely long dS plateau at ip a, 
exponentially rapidly approaching its asymptotic value 


V dS = ap 2 . 


(3.3) 


Predictions from such theories provide a very good fit to 
observational data for a broad class of functions f(Z) as 
discussed in mm- 


However, for such analysis to hold, it is important to verify 
that S = se lJ = 0, and x = Im Z = 0, or to find a way to 
stabilize these fields at their zero values. The point S = x = 
0 is indeed an extremum of the potential for S and x, but 
one should also check whether this extremum is a minimum, 
or a maximum of the potential. 


Let us start with the field x. One can show that its mass 
squared is given by 

™ l ( z ) = !( p ( 6 a / 2 ( z ) + I 1 - z2 ) 2 [(/'(-)) 2 - /( 2 )/"( 2 ) l ) ■ 

( 3 - 4 ) 

If we consider potentials V = ap 2 f 2 (z) vanishing at z = 0, 
then at the minimum one has /(0) = 0, and m 2 ( 0) is positive 
and coincides with the inflaton mass squared at that point, 

m 2 x (0) = m 2 M = ll(fm 2 }) ■ ( 3 - 5 ) 
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Meanwhile (3.4) implies that during inflation 


\{*) = 2 V{z) = 6 H 2 


(3.6) 


up to small corrections suppressed by the slow-roll param¬ 
eters. Thus during inflation this field is strongly stabilized 
for all values of a and tends to rapidly roll down to x = 0 
and stay there. 


The condition 5 = 0 can be satisfied e.g. if the goldstino 
superfield is nilpotent m, or if it is strongly stabilized at 
S = 0 by terms ~ (SS) 2 which can be added to the Kahler 
potential EU, 0. Indeed, one of these two options is com¬ 
pulsory in many theories where the Kahler potential has the 
term SS under the logarithm. However, in the theories with 
the Kahler potential (2.2) considered in this section the situ¬ 
ation is much better: The mass squared of the field S = s e l 7 
at s = x = 0 does not depend on 7 and is given by 


m 2 s (z) = ^(l-z 2 )\f(z)) 2 >0. 


(3.7) 


This mass becomes much smaller than H during inflation 

at z ~ 1 , but the point s = 0 is always the minimum of 

the potential with respect to the field S. This field cannot 

2 

change much because of the overall fast growing factor e s in 
its potential. Thus, the classical field s rolls down towards 
s = 0 and stays there even without adding any stabilization 
terms to the theory. Since its mass is small, it can experience 
small inflationary perturbations in the vicinity of s = 0 , 
but these perturbations are inconsequential unless one makes 
some radical assumptions about strong suppression of the 
probability of decay of the field S in the process of reheating 

m- 


The same models can be discussed in terms of half-plane 
variables, with Z => 777 - Thus, one can use the Kahler 
potential ( 2 . 10 ) and 


V 





FIG. 1. T-model potential V = cx.^ 1 2 tanh 2n for a = at 

Im Z = 0. It is plotted in units afi 2 = 1. The potential has two 
symmetric shoulders. 



FIG. 2. Here we present the cosmological observables n s and 
r for simple T-models with different monomial potentials V = 
a/i 2 ^ tanh((/?/-\/6a)^ with n = ( 1 / 2 , 3 / 4 , 7 / 8,15 3/2,2,3) starting 
from the right, increasing to the left, with the purple line for n = 1 
in the middle. We use the logarithmic scale in r. 


WT eral - T =VSnSf(^). (3.8) 

In what follows we will give some simple examples of these 
models, both in disk variables and in half-plane ones. 


4. T-MODELS 


The simplest a-attractor T-model has Kahler potential 

(4.1) 


(2.2) and superpotential (3.1) with f(Z) = Z\ 
W„ simple - T =Va^S Z . 

The inflationary potential for the canonical field ip is 


V = a/j 2 tanh' 


2 V 


Vfkr 


(4.2) 


We show these potentials in Fig. 1 for few values of a. 


One may also consider a more general case where f[Z] = 
Z n and 

V = a/i 2 tanh 2 " —p= ■ (4-3) 

V6a 


We illustrate the attractor properties of these models in 
Fig. 2 in the logarithmic scale in r. The attractor line, start¬ 
ing at r = 10 -3 towards smaller r, remains universal, the 
value of r depends linearly on a and looses any dependence 
on n: 


1 - n s = 



12a 

ivf 


(4.4) 


where N e is the number of e-foldings. In the simplest T- 
rnodels (4.2) the amplitude of the scalar perturbations also 
does not depend on a and matches the Planck 2015 normal¬ 
ization for /x ft: 10 -5 [61. 


Same models i n half- plane variables are found by using the 
Kahler potential (2.10) and replacing Z => in equations 
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of the previous section. In particular, the model with f(Z) = 
Z becomes 


tt t simple-T 

vv m 


y/a p S 


T - 1 
T+ 1 ' 


(4.5) 


5. E-MODELS 


Another interesting class of models, called E-models, ap¬ 


pears if one uses K® (2.2) and the following superpotentials 


tt rsimple-E _ / 

W B p = yj a pS 


2Z 
Z + l 


(5.1) 


The inflationary potential for the canonical field ip is 


In more general case 

^o general - E = VapS /(^) 


and the potentials are 


V = ap 2 f 2 (1 



(5.3) 


(5.4) 


Let us take a simple example of a general case where 

fim = [m n “d 


V 


= ap 2 


( 


1 — e 



(5.5) 


We illustrate the attractor properties of these models in 
Fig. 4 using a logarithmic scale in r. The attractor line, 
starting at r ~ 5 • 10~ 4 towards smaller r, remains universal, 
the value of r depends on a linearly and becomes completely 
independent of n: r = a 



During inflation, in all of these models, the effective mass 
of the inflaton partner, taking into account the non-minimal 
kinetic term, is given by 6H 2 , up to slow roll parameters, 
where H is a Hubble parameter during inflation. Thus in this 
class of improved a-attractor models in supergravity there 
is no need for stabilization terms. 


The models in half-p lane variables are found by using the 
Kahler potential (2.10) and replacing Z => in equations 
of the previous section. In particular, the model with f(Z) = 
becomes 


< mP ' e - E = yf^pS 


T - 1 
T 


(5.6) 


6. ASYMMETRIC T-MODELS 


FIG. 3. Potentials of simple E-models V = (l — e for a = 

1, 2, 3, 4, 5. It has an inflationary plateau for large positive ip and a wall 
at small negative tp. 



In section [4j we considered the simplest T-models which 
have potentials with two shoulders of equal height. In com¬ 
parison, E-models have only one inflationary plateau, with 
the potential infinitely rising at large negative values of the 
field ip, without forming a second inflationary plateau. In 
this section we will consider an intermediate class of mod¬ 
els, where the potentials have two shoulders with differing 
heights [61 [261127]. In a certain sense, one can think about 
such models as interpolating between the simplest T-models 
and E-models, but they are different because each of these 
two asymmetric shoulders may serve as an independent in¬ 
flationary plateau. Such models offer additional flexibility, 
which, as we will see, may be quite handy. 


Let us use disk Kahler potential (2.2) and a superpotential 


in (3.1) with f[Z\ = (1 — c Z): 

WD inear = y/a n S (1 -cZ) . 


( 6 . 1 ) 


FIG. 4. Here we present the cosmological observables n s and r for E- 

models with different monomial potentials V = (1 — e \^ 3a¥> ) 2n with 
n = (1/2, 3/4, 7/8,1, 3/2, 2, 3) starting from the right, increasing to the 
left, with the line for n = 1 in the middle. The scale in r is logarithmic. 


The resulting potential in the inflaton direction in Z vari¬ 
ables is 


V = ap 2 {l -cZ) 2 . 


( 6 . 2 ) 
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We will assume for definiteness that c > 0. The potential 
has a minimum at Z = 1/c, which belongs to the required 
range \Z\ < 1 for c > 1. In terms of the canonical field ip, 
the potential is 

V = ap 2 tanh — 1 ^ . (6-3) 

At large negative values of the inflaton field, the potential 
has a “shoulder” at V- = (c+1) 2 , and at large positive ip the 
potential has a shoulder of a different height, V + = (c — l) 2 , 
so that 


Y± 


c — 1 

c +1 


2 


(6.4) 


This ratio can be arbitrarily small for 0 < c — 1 -C 1. 

The situation becomes even more interesting in the theory 

W„ xp = VapS( l-e~ &z ) (6.5) 

with /3 » 1. Exponentially suppressed terms e~^ z may 
appear e.g. due to non-perturbative effects [42] • The po¬ 
tential in this case has a minimum at Z = ip = 0. It has 
two shoulders of different height, each of which is capable of 
supporting inflation, with the inflaton potential 

V = ap 2 (l - e -/3 tanh2 n vfa) 2 . ( 6 . 6 ) 


Independently of a, for (3 1, the two shoulders have rela¬ 

tive height V+ /V- ~ e -2 ' 3 . Thus for large (3 the potential has 
two shoulders of exponentially different heights, see Fig. [5]for 
a particular case of a potential with a = 1, (3 = 2. 


V 



FIG. 5. Asymmetric shoulders. 


Despite the exponential difference in the height the values 
of the inflationary parameters n s and r for inflation at each 
of these two shoulders coincide and are given by (4.41, at 
least for not too large values of a. However, the amplitudes 
of scalar perturbations produced at these two branches are 
exponentially different from each other. At the upper branch 
the amplitude of the perturbations is proportional to — 
1 ), and at the lower branch it is proportional to p(l — e~@), 
thus smaller by a factor of e _/3 . 


Suppose that the last stage of inflation, which determines 
the large scale structure of our universe, occurs at the lower 
shoulder with large positive ip. Then, according to [ 6 ], one 
should take 


V + /a w p 2 « 1(T 10 (6.7) 

in Planck units. Considering a = 0(1), one finds V+ ~ 
10 -10 . However, the left shoulder can be arbitrarily high, 
e.g, by taking (3 ~ 11 one can easily have VL = 0(1) in 
Planck units in this model. This fact will be important for 
the discussion of the problem of initial conditions in this 
scenario. But before discussing the cosmological evolution 
in these models, we will return to the simplest symmetric 
T-models and study their potentials more attentively. 


7. HYPERBOLIC GEOMETRY AND PROPERTIES 
OF T-MODEL POTENTIAL 


Consider the simplest T-model (4.1), (4.2), see Fig. 1. We 


would like to describe this model in a more detailed way, 
which should help us to analyze initial conditions for inflation 
and the probability that it will take place in this model and 
its generalizations. In order to do so, we will the inflationary 
potential in a form most suitable for our investigation. 

In terms of Z = z + ix the metric of the moduli space, 
which determines kinetic terms, and the potential are given 

by 


ds 2 = 


3a 


(1 - ZZ) 2 


dZdZ = 3a 


dz 2 + dx 2 
(1 — z 2 — x 2 ) 2 


(7.1) 


r = ( 7.2) 

The existence of the flat inflaton direction is not obvious if 


one is looking at the potential (7.2) in the original variables 


0 and x, see Fig. [ 6 ] That is why one should try to represent 
the potential in terms of more adequate variables. 

One can understand the situation better by making a 
change of variables z = tanh^=. For x = 0, the field ip 
plays a role of a canonically normalized inflaton field, and 
the existence of the inflationary shoulders of the potential 
becomes manifest in the variables ip, x, see Fig. [7} 

By looking at Fig. [7] one could get a wrong impression that 
the potential in the vicinity of the inflationary trajectory at 
x = 0 is incredibly steep: it looks like a gorge which becomes 
more and more narrow at large \ip\. However, this is just an 
illusion. As we have found, the curvature of this potential 
in the direction orthogonal to the inflationary trajectory is 
given by 2 V(ip) = 6 H 2 , which is almost exactly constant 


during inflation, see (3.6). We found this result by taking 


into account that the field x is not canonically normalized. 
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V 


FIG. 6. Inflation, in terms of the variables 2: and x , takes place in 
the two corners of the potential. The existence of the inflationary flat 
directions for \z\ — 1 <C 1 is not apparent in these variables. 



FIG. 8. The T-model potential in terms of the variables <p and $ has 
two infinitely long dS valleys of constant width. 


v 


FIG. 7. A narrowing trench in the potential in terms of tp and x 
variables. We see here the same effect as a decreasing size of angels and 
devils towards the boundary of the Poincare disk in Escher’s picture 
Circle Limit IV. However, this narrowing of the trench is just an illusion, 
which disappears when one plots the potential in proper coordinates, 
as shown in Fig. [8] 



* *0 


This potential is shown in Fig. [8} It has a minimum at d = 
0 where the kinetic terms of both fields become canonical, 
ds 2 —► \{dp 2 + dd 2 ) at d —> 0. At large values of tp where 
tanh -J= approaches 1, the plot of the potential in terms of 
<p and d has a dS valley of constant, (^-independent width, 
instead of the rapidly narrowing gorge shown in Fig. [7] This 
fact will be very important for us shortly, when we will study 
the cosmological evolution of the fields tp and d and initial 
conditions for inflation in these models. 


For a better understanding of the structure of this po¬ 
tential, it is instructive to simplify even a little further the 
superpotential of our simplest T-model: Instead of W = 


fapS Z (4.1), let us consider a ^-independent superpo¬ 
tential 


W = ^fapS. (7.5) 

The potential in this model in the ip and d variables is 


Now we will make this conclusion manifest by plotting the 
potential in terms of more adequate variables. 


We will use the <I> variables as shown in eq. (2.5) with the 
kinetic terms in (2.8) The potential in these variables is 


V = ap 2 


tanh 


ip + id 





(7.3) 


where I tanh 1 2 = tanh . tanh _ We ma y a l so 

I v6o= l yoct v6q! 

present it in the form 


V = ap 2 


cosh 



ip — cos 



cosh 



ip + cos 





-3a 

■ (7-4) 


V = ap 2 



(7.6) 


Note that this potential does not depend on the inflaton field 
ip, and has a dS minimum V = ap 2 at d = 0. It represents 
an infinite (^-independent dS valley as shown in Fig. [9] 


As one can easily check, the shape of this valley coin¬ 
cides with the shape of the dS valley in the simplest T- 
rnodel (7.41 in the large ip limit. This potential is manifestly 
shift-symmetric with respect to the field p. It is singular at 


cos - 


j^d —> 0, but this singularity disappears if one uses 


canonical variables % defined by d\ = 




. In the limit 


cos y ^ <C 1, which corresponds to F > otp 2 , the potential 
of the field d in terms of the canonically normalized held \ 
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FIG. 9. An infinite (/^-independent dS valley ( |7.6[ ) in the theory tLH 
in (/?, $ coordinates. 


FIG. 10. Two-shoulder potential in the theory \6.b\ in (/?, $ coordi¬ 
nates. 


8. OTHER EXAMPLES OF a-ATTRACTORS 


acquires the following simple form: 


V (x) = ap 2 



(7.7) 


For a > 1, this potential grows exponentially fast, and there¬ 
fore the field \ should rapidly fall to the minimum of the po¬ 
tential and stay there. Meanwhile the situation for a <C 1/3 
is quite different. It is well known that the potentials of that 
type are su fficie ntly flat to support inflation |2S]. Indeed, in 
the theory (7.7) the universe expands as t 1 / 3 ". The defin¬ 


ing condition for inflation is H <C H 2 , which means that 
the Hubble parameter (and the size of the horizon) does not 
change much during the Hubble time H^ 1 . This condition 
is satisfied for a <C 1/3. Thus for a <C 1/3 inflation can 
begin at very large values of V(ip, fl) corresponding to large 
values of and then the universe inflates while the field x 
slowly rolls down to the bottom of the dS valley. After that, 
inflation becomes driven by the field ip. 

The existence of the dS valley of a constant, (^-independent 
width will play a crucial role in our discussion of initial con¬ 
ditions for inflation. Therefore we would like to emphasize 
that the existence of such dS valleys is not just a property 


of the simplest T-model (4.1), (7.4), but a general property 


of the class of models (2.2), (3.1) studied in this paper. In 


particular, the potential of the asymmetric T-models (6.3), 


(6.5) exhibits two infinite dS valleys of constant width but 


different height, as shown in Fig. [10] Meanwhile in E-models, 
there is one infinite dS valley at ip > 0. 

In the next sections, we will describe the cosmological 


evolution in the simplest T-model with potential (7.41 for 
a > 1/3, as well as for a <C 1/3. 


Before moving further to the discussion of the cosmologi¬ 
cal evolution in this class of models, we will consider another 
class of a-attractors, with a slightly different Kahler poten- 
), with SS placed under the sign of 

the logarithm: 


tial as compared to (2.6 


K = —3a log 


cosh 


$ - <h\ 
\J6 a ' 



( 8 . 1 ) 


We will consider the simplest superpotential W = 
yW p S Z = sja p tanh . In this case, the scalar potential 
is 



Note 

( |7~4| ). 

(cos 


/ p \1-3 a 

cos \ — 1 ?) . (8.2) 

V 3a / 


that this potential almost exactly coincides with 
The only diffe rence is that instead of the term 

/ I .1 Q 

This 


-3a 


in (7.41 now we have (cos J ^ 


1—3a 


leads to the stability condition a > 1/3. For a > 1/3 the 
potential has a dS valley very similar to the one shown in 
Fig. pH with curvature given by = 6H 2 (1 — ^). For the 
special case a = 1/3, the potential looks like a disk jT5], see 
Fig. [IT] in Section [l2| 


9. INFLATION IN A HOMOGENEOUS UNIVERSE 


In this section we will study the evolution of the fields ip 
and d during inflation in the FRW background with various 
initial conditions. Before analyzing this evolution numeri¬ 
cally, we first study it analytically in simplified cases, which 
will help us to build intuition instead of fully relying on re¬ 
sults of numerical investigation. 
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Consider the evolution of the held p along the inhationary 
dS valley at •& = 0 in the limiting case p/yfot 1, when the 
potential in the ip direction is constant, V,is = ap 2 , with 
exponentially good accuracy. Suppose that V < 1 but the 
initial velocity of the held ip is extremely large, so that its 
kinetic energy is much greater than the potential energy, 
p 2 /2 O V ( is = ap 2 . Note that p 2 ~ 10~ 10 in the models 
describing inhationary spectra with Planck 2015 normaliza¬ 
tion so for a = 0(1) one should have Vds ~ 10“ 10 . How 
far will the held ip travel from its initial value 4>q in a space 
with energy density p 2 / 2 + VdSi if initially the universe was 
in the Planck energy state dominated by kinetic energy of 
the held ip 2 / 2 = 1? 

The equation of motion of the canonical held ip at ■& = 0 
is ip + 3 Hip = 0, which in our case satishes: 

p + V?>pyJ p 2 /2 + V = 0. (9.1) 

This is simple enough so as to allow a full analytical solution, 
but instead we will solve it approximately, which will be 
helpful for understanding the general situation. 

We will start with initial conditions such that p 2 / 2 = 
Nk V, where Nk is some number much greater than 1. If 
initially the kinetic energy of the scalar held was Planckian 
and Vds ~ 10~ 10 , we have Ny. ~ 10 10 . Then, during the hrst 
stage of the held evolution, one can entirely ignore its poten¬ 
tial energy, so that the equation becomes ip = — ^/ 3/2 p 2 , 
with the solution [29J 

\ = h' = hiT (9.2) 

This kinetic energy dominated regime continues from the 
• 2 

initial state with y- = NkVds to the point where it becomes 
approximately equal to Vds ■ During this time, the held p 
moves by Ap = ^ In Nk ~ 9.4 in Planck units, for Nk ~ 
10 10 . 

After that, the kinetic energy of the held becomes sub¬ 
dominant, and it continues moving in accordance with equa¬ 
tion <p = —py/3Vds. One can show that it passes a distance 
slightly smaller than 1 and stops. Thus the total shift of the 
held p due to its initial velocity is about 10 Planck units, 
independently of a, after which all memory about the ini¬ 
tial velocity of the held at the Planckian time completely 
disappears. This conclusion is confirmed by the results of 
numerical calculations. 

In this investigation, we ignored the motion of the held i). 
A combined motion of these two helds for various a is more 
difficult to analyze, but our numerical calculations show that 
the hnal result remains the same: If one starts with general 
initial conditions corresponding to the Planck energy density, 
including kinetic and potential energy of all helds, and follow 
the evolution of these helds at large p where the potential is 
shift-symmetric with respect to p , then the maximal distance 
the held p may travel from its initial position until it enters 
the slow-roll inhationary regime is Ap = In Nk ~ 9.4. 


In the slow-roll inhationary regime, which starts when 
nearly all of the initial kinetic energy of the held p is lost, 
this held p rolls towards the minimum of the potential by the 
distance Sp = yqf In until inflation ends [T{ 0J . Here 
N e is the number of e-foldings during inflation. For example, 
for a = 1 and N e = 60 one has Sp ss 6 . 

Thus we are coming to a rather interesting conclusion: 
Even if the inhationary plateau is extremely low, compared 
to the Planck density, Vds ~ 10~ 10 , initial motion of the 
helds p and $ does not affect the onset of the inhationary 
regime, at this plateau, if inhation begins at a distance Sp > 
10 from the initial value of the held po at the Planck time. 
For the simple T-models with a = 0(1), this means that 
the last 60 e-folds of inhation do not depend on the initial 
velocities and of the helds p and $ if the initial value of the 
held p was \po\ > Ap + Sp > 16. This conclusion will be 
important for us later, when we will discuss a more general 
situation where the universe in the beginning of its evolution 
may be very inhomogeneous. 


To illustrate this conclusion, we present here some of the 
results of our numerical calculations for a particular case 
a = 1 , using a code developed in m- Using the kinetic 
term in (2.8), and the potential in (7.41 given as a function 
of p and $, we solve the Friedmann equations for the time 
evolution of scalars for a set of initial conditions. We plot 
the trajectories on a contour plot of the potential, the values 
of the helds as a function of time and the scale factor as a 
function of time. The examples here involve various choices 
of initial positions, as shown in Fig. [Tl] and vanishing initial 
velocities. 


To develop some intuition, we begin with the case with 
a = 1 where the scalar helds start with zero initial velocities 
at the height 10 times above the bottom of the dS valley with 
Vds = ap 2 . This will show us what the helds ’want’ to do if 
not pushed hard, and whether this may lead to inhationary 
regime. 


Figs. 11 12 13 and 14 show 12 different trajectories with 
different initial values of the helds such that U(</>o, 0 o) = 
10a/i 2 . The initial positions of the helds are shown by small 
circles. In all cases, the held •& rapidly rolls down to the 
bottom of the valley, oscillates and relaxes there. After that, 
for V^o <1 3 (in Planck units), the universe enters the stage of 
inhation driven by the held p. This stage is longer than 60 
e-foldings for p 0 > 6 . 


Examples above and the corresponding Figures had same 
value of the initial energy density due to a proper choice of 
the initial values of po and at t = 0. However, the initial 
values of velocities were chosen to be zero, and the initial 
values of the potential energy was sub-Planckian. Here we 
would like to analyze a situation where the total value of 
the initial energy, kinetic and potential, is of a Planckian 
scale. Namely, here we will start with initial condition with 
velocities and positions such that at a( 0 ) = 1 the initial 
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FIG. 13. The field ip remains nearly constant and very slowly moves 
towards ip = 0 and oscillates there. The cases with smaller initial values 
of ip < 6 reach the minimum earlier and inflate less. 


In a 


FIG. 11. Time evolution of scalars on a contour plot of the potential. 
In the beginning, the field $ moves towards the minimum of the dS val¬ 
ley, whereas the field ip remains nearly constant, if its initial value was 
large. Then after a short stage of oscillations, the field $ vanishes, and 
the cosmological expansion becomes determined only by the inflaton 
field ip, which slowly moves towards ip = 0 and oscillates there. 


■Q 



FIG. 12. The field $ moves quickly towards the minimum, at the 
bottom of the dS valley. Then after a short stage of oscillations, the 
field $ vanishes and remains in all cases at the bottom of the dS valley. 


energy in Planck units M p = 1 is 


E 0 


= El 


Vn = 1, 


(9.3) 



FIG. 14. The log of a scale factor is plotted, as a function of time. The 
models with smallest, in our examples, initial values of ip < 6 produce 
less than 50 e-foldings, however, all other cases with ip > 6 easily lead 
to more than 60. 


where cos y vanishes, the potential (7.4) is extremely 


steep and easily reaches the Planck values V = 0(1). A 
more proper way to understand the growth of the potential 
at large i9 is to replace i9 by its canonical counterpart as 
we did in (7.7). The singular growth of potential close to the 

point where cos = 0 is replaced by exponential growth 

at large %• For a <1 1/3 the potential (7.7) is very steep 
and field \ rapidly falls towards the dS valley. However, for 
a<l/3 the potential is sufficiently flat to allow inflation at 
all sufficiently large values of V V,is ~ 1CT 10 . 


where 


Ekin — 


If 2 + 1/ 


2 cos 2 




(9.4) 


and the potential is given by (7.4). For the numerical calcu¬ 


lations we took the height of the dS valley of the potential 
at d = 0 and large ip as ap 2 = 10~ 10 in Planck units. 

Our potential has a dS valley at a very small height Vds ~ 
IQ -10 in Planck units, however, when •& is close to the value 


We start at Planckian energies En by making random 
choices of •&, tp, $, (p such that Eq. (9.31 is satisfied and solve 


all equations numerically. Our goal is to find the trajectory 
of the fields starting from the Planckian total energy and 
find out how they reach the bottom of the dS valley loosing 
the total energy by 10 orders of magnitude. We show a sam¬ 
ple of such trajectories on the contour plot of the potential 
in Fig. |15[ For illustration purposes, we took all initial con¬ 
ditions at ip > 9 to illustrate what happens at sufficiently 
large p along the infinitely long dS valley, where the poten¬ 
tial practically does not depend on ip. 
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FIG. 15. Evolution of the fields tp and r) starting at the Planck energy 
density.The fields start on the wall with near maximal i) or anywhere 
between the wall and the bottom of the dS valley with velocities such 
that the total energy is always Eft + Vo — 1. In this process their 
total energy is decreasing to ~ 10" 10 when they reach the bottom of 
the valley. 


which just starts to move to smaller values of ip. 


In Fig. 16 we show the time evolution of the Hubble pa¬ 


rameter. For all initial conditions in Fig. 15 the Hubble 
constant H starts at 3Hq = E 0 = 1, i.e. H 0 = ^, and ends 
at the bottom of dS valley with V ~ 10” 10 . The interme¬ 
diate stage is somewhat different for examples with various 
positions and velocities, however, eventually they all end up 
with 3 H 2 « Vds ~ 10 —10 , H « -AlO” 5 . The intermediate 
flattening of the Hubble parameter shown in the Figure has 
an interesting origin: it is due to d reaching its maximum, 
before it goes down to zero. 


As we see, if we start with sufficiently large initial value of 
the field ipo > 10, then, independently of the choice of other 
initial conditions at the Planck density with E% in + Vo = 1, 
the fields fall to the bottom of the valley and slow down 
there in such a way that the field tp does not change much in 
this process, in complete agreement with our general expec¬ 
tations. Thus, all choices of initial conditions with arbitrary 
velocities of the fields and arbitrary deviations of •& from zero 
(from the bottom of the valley) will lead to inflation if the 
initial value of the held ip took any value in the infinitely 
long interval 10 < |^jq| < oo. 


H 
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So far, we studied the cosmological evolution in a homoge¬ 
neous universe, with the same values of the fields everywhere. 
Now we are going to discuss a more general case where the 
universe may be non-uniform. Before doing it, we will give a 
brief review of the problem of initial conditions for inflation 
in an inhomogeneous universe, and later on we will return 
to the cosmological attractors in this context. 


10. INITIAL CONDITIONS FOR INFLATION IN 
AN INHOMOGENEOUS UNIVERSE: A BRIEF 
REVIEW 


FIG. 16. The Hubble parameter for all initial conditions in Fig. |15| 
starts at 3 Rq = 1, Hq = A=, and ends at the bottom of dS valley with 
V ~ 10" l0 . The intermediate stage is somewhat different for examples 
with various positions and velocities. 

The fields start on the wall with near maximal d or any¬ 
where between the wall and the bottom of the dS valley with 
velocities such that the total energy is always E% in + Vo = 1. 
In this process their total energy is decreasing to V,is ~ 10” 10 
when they reach the bottom of the valley. It may not be ap¬ 
parent from this figure, but after a series of oscillations with 
a rapidly decreasing amplitude, similar to what was shown 
in Fig. [T3], the fields relax at the bottom of the valley. For 
large values of p, the potential along the valley is almost ex¬ 
actly flat. Therefore the inflaton field stay almost still for a 
long time until it start slowly rolling to the minimum of the 
potential at tp = 0. The beginning of this slow motion is only 
barely seen here in the end of the lower (yellow) trajectory, 


The theory of initial conditions for inflation was developed 
simultaneously with the invention of the chaotic inflation 
scenario HBIHSH221, see [23] for a recent discussion. In the 
previous sections we discussed only the homogeneous uni¬ 
verse case, where the situation is relatively simple. Now we 
are going to consider the initial condition problem in a more 
general situation, when the universe may be inhomogeneous. 

Let us remember first what was the main problem with 
the hot Big Bang scenario in this respect. In that model, 
the universe was born at the cosmological singularity, but it 
was possible to describe it in terms of classical space-time 
only when time is greater than the Planck time t p ~ 1 . 
At that time, the temperature of matter was given by the 
Planck temperature T p ~ 1, and the density of the universe 
was given by the Planck density p p ~ 1. The size of the 
causally connected part of the universe at the Planck time 
was ct p ~ 1 . Each such part contained a single particle with 
the Planck temperature. The subsequent evolution was sup- 
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posed to be adiabatic. As such, the total number of such 
particles is approximately conserved, so it must be greater 
than the total number of particles in the observable part of 
the universe, n ~ 10 90 . This means that the universe at the 
Planck time consisted of 10 9 ° causally disconnected parts. 
The probability that all of these independent parts emerged 
from singularity at the same time with the same energy den¬ 
sity and pressure is smaller than e ~ 10 , and the probability 
that our universe from the very beginning was uniform with 
an accuracy better than 10 -4 is even much smaller. 

The change of the perspective brought by the invention 
of the chaotic inflation scenario was quite dramatic. The 
main condition required in the simplest models of chaotic 
inflation was the existence of a single Planck size domain 
where the kinetic and gradient energy of the scalar held is 
few times smaller than its potential energy V(<j>) ~ 1. For 
sufficiently hat potentials, it leads to inhation, so the whole 
universe appears as a result of expansion of a single Planck 
size domain. According to naima, the probability of 
this process is not exponentially suppressed. After that, the 
universe enters an inhnite process of self-reproduction m- 

However, in the argument given above it was assumed 
that inhation may start at the densities comparable with 
the Planck density. This condition is not satished in hilltop 
models [3T], in the Starobinsky model j7], in the Higgs in¬ 
hation model uni, and in the simplest versions of the broad 
class of the cosmological attractor models discussed in this 
paper and in HE], where inhation occurs only at V <C 1. 
Thus the issue of initial conditions in these models should 
be addressed. 

To explain the main problem, we will consider the new in- 
hation/hilltop models [31] and follow the argument against 
them and in favor of chaotic inhation given in [22] . Inhation 
in such models begins at the density about 10 orders of mag¬ 
nitude below the Planck density. If the universe is closed, 
then a typical closed universe born at the Planck density and 
filled with relativistic particles collapses within the Planck 
time, unless it is extremely large from the very beginning 
[221 . But then inhation does not explain why our universe is 
large. Rather, inhation in such models requires the universe 
to be large. 

Moreover, it is hard to expect a large universe to be uni¬ 
form on the scale 7L _1 , as required for inhation, at the time 
when the energy density of normal matter drops down to 
the energy density at the top of the inhationary potential 
V ~ 10 -1 °. Indeed, the size of the Hubble domain at that 
time is about t ~ 10 5 , but the total size of the hot uni¬ 
verse was growing only as y/t. As a result, the Hubble size 
domain at the moment when inhation could start contains 
t 3 / 2 ~ 10 7 parts which had the Planck size and were causally 
disconnected at the Planck time. The probability that the 
universe in all of such domains had the same density can 
be estimated by e ~ 10 . It is much better than e -10 , but 
it is still highly unsatisfactory [221 - Recently, a very simi¬ 


lar argument against naturalness of initial conditions in the 
Starobinsky model was given in [32] . Does it mean that all 
such models are in trouble? 

Fortunately, several ways to solve the problem of initial 
conditions in the models with low scale of inhation have been 
proposed more than 10 years ago, see a discussion of this 
issue in [25] and references therein. One of the simplest 
solutions is to consider a hat compact universe having the 
topology of a torus, Sf, 

ds 2 = dt 2 — a 2 (t) dx 2 (10.1) 

with identification Xi + 1 = Xi for each of the three dimen¬ 
sions [33] . Whereas this may sound exotic, such models are 
routinely used by string theorists who assume that 6 inter¬ 
nal dimensions are compactihed, and by experts in numerical 
relativity who study the universe on a lattice with periodic 
boundary conditions. The trick is to make this box expo¬ 
nentially big and its periodicity unobservable, which can be 
easily done by inhation. 

Suppose, for simplicity, that a± = <12 = <13 = a(t). In this 
case the curvature of the universe and the Einstein equa¬ 
tions in terms of a(t) will be the same as in an inhnite hat 
Friedmann universe with the metric ds 2 = dt 2 — a 2 (t)dx 2 . 
In our notation, the scale factor a(t) is equal to the size of 
the universe in Planck units. 

Let us assume that at the Planck time t p ~ 1 the universe 
was radiation dominated, h < T 4 = 0 ( 1 ) and that at the 
Planck time the total size of the box was Planckian, a(t p ) = 
0(1). Then the whole universe initially contained only 0(1) 
relativistic particles such as photons or gravitons, so that the 
total entropy of the whole universe was 0 ( 1 ). 

The size of the universe dominated by relativistic parti¬ 
cles grows as a{t) ~ y/t, whereas the mean free path of the 
gravitons grows as H _1 ~ t. If the initial size of the universe 
was 0 ( 1 ), then at the time 1 > 1 each particle (or gravita¬ 
tional perturbation of the metric) within one cosmological 
time would run all over the torus many times, appearing in 
all of its parts with nearly equal probability. This effect, 
called “chaotic mixing,” should lead to a rapid homogeniza¬ 
tion of the universe [551KT5] . 

In fact, to achieve a modest degree of homogeneity we do 
not even need chaotic mixing. Indeed, density perturbations 
do not grow in a universe dominated by ultra-relativistic par¬ 
ticles if the size of the universe is smaller than H~ x . This 
is exactly what happens in our model before inflation be¬ 
gins. Therefore the universe should remain relatively homo¬ 
geneous until the thermal energy density drops below Vds- 
Once it happens, the universe rapidly becomes exponentially 
large and homogeneous due to inflation. 

The only issue with this scenario is that at the time when 
the thermal energy density drops below Vds , the scalar field 
should be at the position where the potential V(<j>) is flat and 
inflation is possible. This is not easy to achieve in the new 
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inflation/hilltop inflation scenario because the probability to 
land by chance at the tiny vicinity of a maximum of the 
potential in these models is not very high. However, in the 
models with long inflationary slopes of the potential, such 
as the simplest chaotic inflation models m or especially in 
the models discussed in this paper, where the inflationary 
slope of the potential is infinitely long, this last part should 
be very easy. 

Thus in this scenario, just as in the simplest chaotic infla¬ 
tion scenario, inflation may begin if we had a sufficiently ho¬ 
mogeneous domain of a smallest possible size (Planck scale), 
with the smallest possible mass (Planck mass), and with the 
total entropy 0(1). We see no reason to expect that the 
probability of formation of such universes is strongly sup¬ 
pressed. 

One can arrive at a similar conclusion from a completely 
different point of view. Investigation of the quantum cre¬ 
ation of a closed or an infinite open inflationary universe 
with b <C 1 shows that this process is forbidden at the clas¬ 
sical level, and therefore only occurs by tunneling. The main 
reason for this result is that closed de Sitter space always 
has its size greater than H~ x ~ 1/y/V, and the total energy 
greater than H~ 3 V ~ 1/y/V. Than is why the universe with 
V <C 1 is large, heavy and difficult to create. As a result, 
the probability of this process is exponentially suppressed 
mm- Meanwhile, creation of flat or open compact uni¬ 
verses is possible without any need for tunneling, and there 
would be no exponential suppression for the probability of 
quantum creation [33 HE H7]. This suggests that the prob¬ 
lem of initial conditions for low energy scale inflation can be 
easily solved if one considers topologically nontrivial com¬ 
pact universes. This possibility is completely consistent with 
observations, since inflation lasting more than 60 e-foldings 
makes all topological effects unobservable. 

Thus we believe that there is no real problem with ini¬ 
tial conditions for inflation in a broad class of inflationary 
models even if inflation may occur only at V <C 1. Other 
ways to solve the problem of initial conditions for low en¬ 
ergy scale inflation, including the possibilities related to the 
string landscape scenario pMFHj . can be found in [23j. As 
we will see now, the cosmological attractor models studied 
in this paper have additional advantages in this respect. 


11. COSMOLOGICAL ATTRACTORS, DS VALLEY, 
AND INITIAL CONDITIONS FOR INFLATION 


In order to make our basic idea clear, we would like to 
start with a discussion of a simple toy model of a massless 
scalar field with gravity and a cosmological constant V = 
const -C 1: 


(n.i) 


The description of the evolution of the universe in terms of 
classical space-time begins when the sum of the kinetic and 
gradient energy of the field <p and the vacuum energy V for 
the first time becomes equal to the Planck energy density: 
\{p 2 + ditp 2 ) + Vds = 1- For Vds 1) this means that 
§(^ + 0 i¥ > 2 )«l. 

What happens with such a universe in the future? If it 
is a compact open or flat universe, then we believe that 
the answer is already given in the previous section. The 
field <fi is massless, its perturbations represent relativistic 
waves/particles, and because of their ultra-relativistic equa¬ 
tion of state combined with the chaotic mixing, the uni¬ 
verse stays small and relatively homogeneous, of the size 
a(t) -C H -1 . Its size grows as y/t until the energy density 
of the scalar field p v (t) ~ t~ 2 drops below V. After that, 
the universe becomes an exponentially expanding dS space, 
which will contain perturbations of the scalar field with an 
exponentially decreasing amplitude. No other choice is avail¬ 
able for the system. 


One may consider a more complicated case: the same 
model (11.1), but the universe was born large from the very 
beginning. In our opinion, this assumption is problematic 
and unnecessary. If everything that we see now could origi¬ 
nate from something small, why should we consider a much 
more complicated way to reach the same goal? One of the 
problems addressed by inflation was to explain its huge size. 
Since we know how to solve this problem using inflation, an 
assumption that the universe was large prior to inflation is 
no longer required. Nevertheless, for the sake of argument, 
and to cover all options, we will consider this possibility as 
well. Then there are two basic possibilities: 


1) The whole universe collapses before the onset of ex¬ 
ponential expansion in dS space with the cosmological con¬ 
stant V. But this is possible only if the collapse occurs very 
quickly, before the energy density of the initially expanding 
universe becomes smaller than Vds and dS inflation takes 
over. For Vds ~ 10~ 10 , this may happen only if the whole 
universe collapses almost instantly, within about 10 -28 sec¬ 
onds from the birth of the universe, which means that we are 
not talking about a universe with human observers trying to 
explain its homogeneity, but about a short-living quantum 
fluctuation. 


2) If we are not talking about a small instantly collapsing 
“virtual” universe, then the next possibility to consider is a 
large expanding universe with initial inhomogeneities, which 
may grow and become large in some parts of the universe, 
but not everywhere. For example, the universe may contain 
such local inhomogeneities as primordial black holes. But if 
the universe as a whole continues to expand, then density 
of matter in all of its expanding parts rapidly drops down 
to Vds- These parts will continue to grow exponentially, 
whereas the collapsing parts of the universe will not. As a 
result, the volume of the universe will become dominated by 
the exponentially large volume of dS space. Even if it con- 
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tains black holes or other large inhomogeneous stable struc¬ 
tures, they will move away from each other exponentially fast 
and therefore eventually they will become irrelevant. This is 
the similar to the fate that awaits all of us if the dark energy 
is a cosmological constant: Eventually all distant galaxies 
and other inhomogeneities will move exponentially far away 
from us, and the universe will become an empty dS space. 

Thus we expect that in the presence of a sufficiently large 
cosmological constant V,is- at least some parts of a large 
expanding universe do not collapse. The collapsed parts are 
of little interest to us, whereas all parts that did not collapse 
within the first 10 -28 seconds join the exponential expansion 
driven by the positive cosmological constant Vas¬ 
in other words, in the context of this simple model, the 
problem of initial conditions all but disappears. Instead of 
asking whether it is probable that a large inhomogeneous ex¬ 
panding universe enters a stage of inflation, one may wonder 
what kind of miracle may prevent at least some parts of the 
universe to enter the process of exponential expansion in an 
eternally growing dS space? 


To explain the relation between this toy model of dS space 
and the theory of initial conditions for inflation, let us take 
the next step and consider yet another toy model (7.51 with 
an exactly shift-symmetric dS valley potential (7.6) shown 
in Fig. [9j This potential represents an infinite dS valley with 
the shape coinciding with the shape of the dS valley in the 
simplest T-model (7.41 in the large p limit. According to 
the argument given above, the universe in the model (7.5) 
should end up at the bottom of dS valley and enter the state 
of exponential expansion in dS state at i9 = 0. 


Note that even after rolling to the dS valley at i9 = 0, 
different parts of the universe may contain different values 
of the field p, depending on the initial conditions at the 
Planck time. Indeed, the theory (7.51, (7.6) is exactly shift- 
symmetric with respect to the Held p. This means that all 
initial conditions for the field p are equally probable. There¬ 
fore all values of this field after it falls to the dS valley are also 
equally probable: A large dS universe becomes populated 
with all possible values of the field (p in its different parts, 
but the gradients of this field rapidly become exponentially 
small due to dS inflation, and the universe becomes locally 
uniform, up to small inflationary quantum fluctuations of 
the field p. 


This is not a complete inflationary scenario yet, because in 
this model the universe is forever trapped in dS space. In or¬ 
der to implement the standard slow-roll inflation, one should 
return to realistic models with Z-dependent superpotentials 
discussed in our paper, such as W = y/a pS Z proposed in 
eq. (4.1). These models were studied in Section [7] In such 
models, the field does not stay in dS valley but slowly rolls to 
the minimum of the potential at small ip. But from the point 
of view of the theory of initial conditions, these models do 
not differ much from the toy models of dS space considered 


above. 


Indeed, the kinetic terms in all of these models coincide 


with the kinetic terms of the theory (7.51, and the potential 


of the simplest T-models (7.4) coincides with the dS val¬ 


ley potential (7.6) with exponentially good accuracy every- 

< 


where except a finite interval of width \ip\ < yfa. According 
to our results, even if the kinetic energy of the field ip ini¬ 
tially was as large as the Planckian energy, this field does 
not move by more than |A(/j| = 10 until it slows down and 
reaches the slow-roll inflationary regime along the dS valley. 
Therefore for an infinitely large range of initial conditions 
v'a + 0(10) < |<£o| < oo, the first stages of the cosmological 
evolution in T-models coincide with the evolution in the sim¬ 


ple dS model (7.5), (7.6) with exponentially good accuracy. 


More generally, one can show that the inflationary regime in 
this class of models occurs if, after the stage of rolling, the 
value of the inflaton field at the bottom of the dS valley re¬ 
laxes at \ip\ > min[0( v / a), 0(1)]. This means that the stage 
of slow-roll inflation occurs for all values of po at the Planck 
density such that 0(10) < |<£o| < oo. 

The phase volume of initial conditions with \p\ < yW + 
0(10) is finite. Meanwhile the phase volume of all possi¬ 
ble inflationary initial conditions in these models is much 
greater: it is infinitely large, with all values with \<p\ 0> 
y/a + 0( 10) being equally probable because of the shift sym¬ 
metry. For all values of the field p in this infinitely large 
interval, the probability of the slow-roll inflation in the cos¬ 
mological attractors studied in this paper coincides with the 
probability that the universe ends up in an exponentially 


expanding dS space in the dS valley model (7.5). In accor¬ 


dance with the arguments given in this section, this means 
that initial conditions for inflation in the family of cosmo¬ 
logical attractors studied in this paper are quite natural. 


12. FURTHER GENERALIZATIONS 


Before concluding our general investigation of initial con¬ 
ditions for inflation for cosmological attractors, we would 
like to point out some special cases where this investigation 
is even simpler, and also to extend our results to other infla¬ 
tionary models. 


The first case involves asymmetric T-models (6.3), (6.5) 
shown in Figs. [5] and [lO] In these models, for appropriate 
choice of parameters, the hight of the upper shoulder (i.e. of 
the dS valley at a higher altitude) can be Planckian. Thus 
inflation may start at the Planck density, as in the simplest 
models of chaotic inflation, in which case the problem of ini¬ 
tial conditions is easily solved [21 IT8U22) . When inflation 
ends at the upper shoulder, the field p rolls down to p = 0 
and, by inertia, climbs the lower shoulder. As we have al¬ 
ready shown, the field cannot move by more than Ap ~ 10. 
After that, the field stops, and the second (last) stage of in- 
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flation begins, which is responsible for the structure forma¬ 
tion in the observable part of the universe. We checked that 
for a = 0(1) the duration of the second stage of inflation is 
more than sufficient to account for the last 60 e-foldings of 
inflation. 


Another interesting possibility emerges in supergravity 


models with the Kahler potential (2.2) and a -C 1/3, as 


well as in the models with the Kahler potential (8.1) and 


a = 1/3 + 6 where 0 < S -C 1/3. In such models, inflation 
can begin directly at the Planck boundary with large val¬ 


ues of the canonical field x i n Eq. (7.7). Then the problem 


of initial conditions for such models is resolved in exactly 
the same way as in other versions of chaotic inflation where 
inflation may start at the Planckian values of the potential 
V = 0(1 ) HS1HSH22]- 

The arguments given above can be easily extended to the 
models discussed in Section[8]with the Kahler potential (8.1). 
As we already mentioned, the potential in these models for 
a > 1/3 also has an infinite dS valley ( |8.2[ ) with the ip- 
independent curvature m| = 6H 2 (1 — ), so most of our 

arguments apply there for the case a > 1/3. The special case 
a = 1/3 is particularly interesting from the point of view of 
the theory of initial conditions, as already discussed in m- 
Indeed, in this case it is convenient to use a different set of 
coordinates Z = e lQ tanh ^ and find that, for the simplest 
superpotential W = y/a fi S Z, the inflaton potential does 
not depend on the angular variable 6 and is given by V = 
afi 2 tanh 2 ^=. 



FIG. 17. The potential V in the theory \8.l\ with W = y/a p, S Z for 
ol = 1/3 in terms of the inflaton field tp and the angular variable 6. 
Everywhere except for the region with ip < 4, the potential is flat with 
exponentially good accuracy, which provides perfect initial conditions 
for inflation m- 

The shape of the potential in these variables is shown in 
Fig. [lT] The last 60 e-foldings of inflation in this model 
are determined by a tiny part of this figure at ip < 4. The 
potential V at large p is almost exactly constant, and the 
phase volume for initial conditions corresponding to p > 4 
is infinitely large. This provides perfect initial conditions for 
inflation. Independently of the initial velocity of the scalar 
fields, their kinetic energy rapidly dissipates due to the cos¬ 
mological evolution. The fields freeze at some point of the 
infinitely large plateau, until the exponentially slow descent 
in the radial direction towards the minimum of the poten¬ 


tial shown in Fig. [lT] begins [T5] . All of the arguments given 
above concerning initial conditions for exponential expansion 
in dS space are immediately applicable to this model. 


Up to now, we concentrated on a special class of a- 
attractors in supergravity, describing a combined cosmolog¬ 
ical evolution of two fields, p and d. We found that in some 
exceptional cases the field d can help inflation to begin (e.g. 
in the model (8.1) with a <C 1/3), but in general it does 
not play much role in the discussion of initial conditions for 
inflation: If the field p initially was very large, then it falls 
to the dS valley at very large values of p independently of 
the behavior of the field d. Our main arguments concern¬ 
ing naturalness of initial conditions for inflation in Section 
EH take this fact into account. Therefore they equally well 
apply to all cosmological attractor models based on a single 
inflaton field with a plateau dS potential, such as (1.1). This 
includes the Starobinsky model [7j, Higgs inflation [lOj, con¬ 
formal attractors [T|, single-field supergravity attractors [3], 
and universal attractors with non-minimal coupling to grav¬ 
ity [2] OS] • More generally, these arguments should apply to 
any model where the inflaton potential has a sufficiently long 
and flat slope, as in the simplest models of chaotic inflation 

Ga. 


13. SUMMARY 


In this paper we described some simple versions of the 
theory of the cosmological attractor models mil Ena Hz]. 
We have used a novel set of coordinates describing the ge¬ 
ometry of scalar Helds. It corresponds to a choice of the 
Killing adapted geometry where the inflaton shift symme¬ 
try, p —> p + const of the Kahler potential, is manifest and 
is only broken by the superpotential. As the result, the po¬ 
tential of these refined models has an interesting universal 
property: at large values of the inflaton field p it looks like 
an infinite dS valley of nearly exactly constant depth and 
width, see Fig. [8] In other words, at large p the poten¬ 
tial has shift symmetry with respect to the shift of the field 
p. This property of the theory is not immediately appar¬ 
ent from its formulation, but can be established using the 
properties of the hyperbolic geometry of the moduli space 
discussed in DUE] and further studied in our paper. The 
novel element contained in this paper is the choice of con¬ 
formally flat coordinates, and also a particular choice of the 
Kahler potential, which makes the shift symmetry manifest 
and ensures stability of the dS extremum. 


We believe that this interesting property is sufficient for 
solving the problem of initial conditions in this class of infla¬ 
tionary models. This can be done in several different ways 
discussed in [23] and in Section 10 of our paper, but the ex¬ 
istence of the shift symmetry revealed in this paper allows 
for yet another solution which is rather simple and general. 
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Our investigation shows that during the pre-inflationary 
stage of expansion of the universe, the canonically normal¬ 
ized inflaton field ip cannot change by more than 0(10). 
Therefore for all initial conditions for the field ip, except 
for a finite interval near the minimum of the potential at 
p = 0, the fields roll down to the lower part of the dS valley 
far from ip = 0 , and inflation almost inevitably begins there. 
This argument is directly applicable as well to all cosmolog¬ 
ical attractor models, such as <0>> based on a single real 
inflaton field (f> with a plateau potential, and to any other 
model where the inflaton potential has a sufficiently long 
and flat slope. 
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